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We present theoretical calculations of the thermal conductivity for the "ac- 
cidental degeneracy" and "enlarged symmetry group" models that have been 
proposed to explain the phase diagram of UPt$. The order parameters for 
these models possess point nodes or cross nodes, reflecting the broken symme- 
tries of the ground state. These broken symmetries lead to robust predictions 
for the ratio of the low-temperature thermal conductivity for heat flow along 
the c axis and in the basal plane. The anisotropy of the heat current response 
at low temperatures is determined by the phase space for scattering by im- 
purities. The measured anisotropy ratio, n c /Kb, provides a strong constraint 
of theoretical models for the ground state order parameter. The accidental 
degeneracy and enlarged symmetry group models based on no spin-orbit cou- 
pling do not account for the thermal conductivity of UPt^. The models for 
the order parameter that fit the experimental data for the c and b directions 
of the heat current are the 2D E\ g and E% u models, for which the order pa- 
rameters possess line nodes in the ab-plane and point nodes along the c axis, 
and the A\ g ®Ei g model of Zhitomirsky and Ueda. This model spontaneously 
breaks rotational symmetry in the ab-plane below T c i and predicts a large 
anisotropy for the ab-plane heat current. 



2 M. J. Graf, S.-K. Yip and J. A. Sauls 

The observations of multiple superconducting phasesSi'i'B of UPt3 have 
led to several Ginzburg-Landau models based on different symmetry groups 
or symmetry breaking scenarios in order to account for the phase diagram. In 
this paper we show that these models for the phase diagram, which are based 
on different pairing symmetries, lead to qualitatively different predictions for 
the thermal conductivity at very low temperatures. 

The pairing models of UPt3 which account for the main features of 
the phase diagram may be grouped into three classes. One class of models 
is based on a two-dimensional representation of Dq^ with the multicom- 
ponent superconducting order parameter coupled to a .Symmetry- Sreaking- 
iueld (SBF). The leading candidates for the pairing state are the even-parity 
Ei g model and the odd-parity model. Both models require a symmetry 
breaking field in order to split the transition in zero field and produce mul- 
tiple superconducting phases in a field. The SBF is generally assumed to 
be the in-plane AFM order parameter that onsets at Tm — 5Kp however, 
other explanations of the SBF have been suggested .BSh 

Several authors have argued that the apparent tetracritical point in 
the H-T phase diagram for field orient ation s gut of the basal plane is in- 
compatible with the 2D E-rep models.li3llJ'll3'll3 This led to explanations 
of the phase diagram in terms of multiple order parameters that are unre- 
lated by symmetry_fe 1 o. '^li?-models'),li3'EJ as well as further examination of 
the E-rep models. ot3 The four possible E-rep models are not equivalent; 
weak hexag onal anisotropy, as is reflected by the weak in-plane anisotropy 
of H C 2 ,03113 allows for an apparent tetracritical point for all field orientations 
provided the order parameter belongs to an E2 orbital representation.EJ The 
odd-parity, E2 U representation with strong spin-orbit locking of the order pa- 
rameter with d| |c (d is the quantization axis along which the pairs have zero 
spin projection, i.e. d-S = 0) also accounts for the anisotropic paramagnetic 
limiting of H C 2 observed at low temperatures. EE The spin-singlet Ei 5 model 
appears to be incompatible with both the tetracritical point for H JL c and 
the anisotropic Pauli limit for H C 2. However, Park and JoyntS argue that 
there is enough freedom in the Ei ff model to account for all existing experi- 
mental data. Analysis of current heat transport data at low temperatures is 
accounted for equally well by either E-rep. We suggest below that ultra-low 
temperature measurements of the thermal conductivity for varying impurity 
concentrations can also differentiate between the two E-rep models. Both 
E-rep models have recently been challenged by the observation of a nearly 
temperature independent Knight shift for H| |cx! Tou and co-workerJll favor 
the ESG-model of Machida et aZ.ta based on spin-triplet pairing with weak 
spin-orbit coupling. Thus, our current understanding of the spin structure 
and parity of the order parameter for UPt3 is unclear because there is con- 
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flicting information from the NMR measurements of the spin susceptibility 
and measurements of the anisotropic paramagnetic limit of the upper critical 
field. However, as we show below, the low temperature thermal conductiv- 
ity of UPt3 is in strong conflict with the non-unitary spin-triplet model of 
Machida et alM 

The j4i?-models are based on an assumed near accidental degeneracy 
of two different pairing channels which are unrelated by symmetry.0'0 In 
these models the phase diagram is accounted for by two primary order pa- 
rameters belonging to different irreducible representations of Dq^, which are 
selected in order to enforce a tetracritical point in the GL theory for the 
H-T phase diagram. The accidental degeneracy models which have been 
investigated in most detail assume that the two order parameters belong to 
different one-dimensional representations of the same parity, e.g. A 2k and 
B lu M Another group of models for the phase diagram may be described 
as Enlarged-Symmetry- Group models. These are hybrids of the accidental 
degeneracy and SBF models. The ESG- models assume an accidental de- 
generacy in the form of a "hidden symmetry", which is lifted by a weak 
SBF coupling, similar to tha t of the E-rep models. Two different ESG- 
models have been proposedJ13E3Ell The model of Zhitomirsky and UedaEil is 
based on an enlarged orbital symmetry group for UPt3. In this model the 
primary pairing channel is assumed to be the d-wave [t = 2) channel of the 
full rotation group, and the degeneracy of the d-wave manifold is lifted by 
weak hexagonal anisotropy: r2[50(3)] — > Ai g © E\ g © Ei g . The Ai g and Ei s 
channels are assumed to be weakly split, in order to account for the double 
transition. Another formulation of an ESG-model assumes spin-triplet pair- 
ing with independent spin and orbital rotation symmetry, i.e. no spin-orbit 
coupling!!! In the model by Machida et a/.0 a one-dimensional, odd-parity 
orbital representation is chosen in order to enforce a tetracritical point for 
all field orientations. The multicomponent order parameter corresponds to 
the three spin-triplet amplitudes described by the d vector, and the SBF is 
assumed to be the AFM order parameter. The main features of each of 
these models of the phase diagram and thermodynamic quantities of UPt3 



are discussed in several articles.1 



In this paper we investigate the low-temperature transport properties 
for the ground states of the AB and ESG models. The heat transport 
coefficients are sensitive to the nodal structure and symmetry of the order 
parameter; at low temperatures the lower-dimensional regions of the Fermi 
surface associated with the nodal regions dominate the electronic heat trans- 
port. The nodes of A(p*) reflect the broken symmetries of the supercon- 
ducting phase. The most specific information regarding the symmetry of the 
order parameter is reflected in the ultra-low temperature limit of the compo- 
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nents of the thermal conductivity tensor, i.e. lim-p^o k/T. The connection 
between the symmetry of the pairing state and the low temperature ther- 
mal conductivity is via impurity-induced Andreev scattering. Impurities 
in unconventional superconductors lead to a band of low-energy Fermion 
states which are formed by impurity-induced Andreev scattering and the 
momentum-space structure of the order parameter £3 In a dilute alloy a 
distribution of impurities leads to a finite density of states, which we call 
an Andreev band. The bandwidth and density of states can be calculated 
from the self-consistent i-matrix and the leading order impurity self-energy 
terms .e3'C3'E3 Impurity scattering within the Andreev band is then the dom- 
inant heat transport mechanism at very low temperatures, i.e. k^T < 7, 
where 7 is the width of the Andreev band. The limiting values for k/T 
are sensitive to the symmetry of the order parameter, the geometry and 
dimensionality of the nodes B^u 

In an earlier papera we compared the electronic thermal conductivity 
measurements on UPt3 by Lussier et aZ.E] with our calculations based on 
the E-rep models. The experimental fact that k c (T)/k c (T c ) < k&(T) / Kb(T c ) 
for temperatures down to T ~ 50 mK provides a strong constraint on mod- 
els for the ground state order parameter and low-energy excitation spec- 
trum. We found excellent agreement between theory and experiment for 
both the Ex 9 and F,2u representation [see also Norman and HirschfeldE3 
and Lussier et aZ.ES] . The quality of our fits for both E-rep models is 
evident in the temperature dependence of the normalized anisotropy ratio 
[K c (T)/K b (T)]/(K c (T c )/K b (T c )] (Figs. 2and3ofRef. |). From the published 
measurements, which extend down to T c /10, it is not possible to distinguish 
between the low temperature phase of an even parity, spin singlet state (Ei 9 ) 
and an odd parity, spin triplet pairing state (E2«). Further experimental 
studies of the dependence of the anisotropy ratio on the impurity concentra- 
tion are needed to distinguish between these two E-rep models. The differ- 
ences in symmetry between the Ej g and E2« models may be tested at ultra- 
low temperatures on nearly clean materials. One possibility for distinguish- 
ing between the Ei s and E2« ground states using transport measurements 
is based on the impurity scattering dependence of the asymptotic values of 
k/T at temperatures below the impurity bandwidth, ksT <C 7. The predic- 
tions for the anisotropy ratio 1Z = ]imT->o[n c (T) / K(,(T)]/[k c (T c ) / K b (T c )] as 
a function of the normal-state impurity scattering rate are shown in Fig. 1 
for the two E-rep models. The key feature is the universal (non-universal) 
limit for the E2 n (Ei s ) model. 
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Fig. 1. The asymptotic low temperature anisotropy ratio of the thermal conductiv- 
ities for an Ei s and E2 M order parameter with the 'best fit' parameters specified in 
Ref. (||) . The arrow indicates the estimated scattering rate obtained from the data 
of Ref. (^ ) for the fit to the E 2u model. 



Below we present our analysis of the electronic thermal conductivity for 
the representative ground states of the AB and ESG-models of UPt3. We 
discuss the variational forms for the order parameter of these models, and 
derive asymptotic limits for the low temperature thermal conductivity. We 
compare our analytical results for the asymptotic limits of k/T with numer- 
ical calculations of the thermal conductivity over the full temperature range 
below T c , and we compare the thermal conductivity calculated within the 
AB and ESG models with the experimental data on UPt3.B^l None of the 
AB or i?5G-models based on spin-triplet pairing without spin-orbit coupling 
account for the temperature dependence and anisotropy of the thermal con- 
ductivity at low temperatures. The only models which can account for the 
thermal conductivity data are the 2D Ei g and F>2 U models, and the orbital 
ESG-model of Zhitomirsky and UedaS The ^4i?-model predicts a large aft- 
plane anisotropy, which differentiates it from the E-rep models. We discuss 
the conditions under which the ^4i£-model can account for the heat transport 
data in UPt3. 



Thermal Conductivity at Low Temperature 

Our analysis focuses on the low temperature region, T < 0.5 T c , where 
heat transport is limited by scattering of electronic quasiparticles off impuri- 
ties. We consider s-wave scattering from a random distribution of defects in 
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the strong scattering limit (5q — ► tt/2). The assumption of nearly resonant 
scattering is in agreement with transport measurements in the normal Fermi- 
liquid and superconducting phases.E^ Inelastic "electron-electron" scattering 
is included by extrapolating the experimentally determined normal-state in- 
elastic scattering rate, l/r ee = AT 2 , to temperatures below T c ; inelastic 
scattering is negligible below T ~ 0.3T C . 

The transport properties of unconventional superconductors are more 
strongly influenced by scattering from non-magnetic impurities than in con- 
ventional superconductors. One of the more striking effects is the ap- 
pearanceof a band of low energy excitations deep in the superconducting 
state.il' S'U'El'ii This occurs for non-magnetic impurities in unconventional 
superconductors when the order parameter changes sign around the Fermi 
surface. Impurity states develop from the combined effects of impurity scat- 
tering and Andreev scattering. In an unconventional superconductor with 
(A(pj)) = and line or cross nodes, a finite concentration of impurities 
leads to Andreev states with a bandwidth, 7 <C Ap. below which the density 
of states is non-zero and almost constant for e < This novel metallic 

band, deep in the superconducting phase, can exhibit universal values for 
the transport coefficients, i.e. independent of the scattering rate at very low 
temperatures, k B T < 7.ii0H§! 

In Ref. pl| it was shown that the principal components of the thermal 
conductivity tensor can be expressed in terms of an effective transport scat- 
tering time that incorporates particle-hole coherence of the superconducting 
state. For T -> 0, 

hm^(T) = V JA ls Trf f (i = a,b,c), (1) 

where v 2 { is the Fermi-surface average of [i>/,j(p/)] 2 , 7s 1 = ^ir 2 k 2 3 Nf is the 
normal-state Sommerfeld coefficient, and 



T 



e//_ f MP/)]V 



^ dpf L"J.n*v^ / ( 2 ) 
2vh J VS rA<W\2-i_,v21§ 



v f,i J [A( P/ ) 2 + 7 

is the effective transport scattering time. The bandwidth of the Andreev 
states is given by 

( 7 (A( P/ )^)-V>) 

cot 2 5 + <7(A( P/ ) 2 + 7 2 )- 1/2 ) 

with the Fermi surface average (. . .), and T u = rii/^nNf). For a given im- 
purity concentration, m, this bandwidth is largest in the limit of unitarity 
scattering, and depends on the order parameter symmetry through the struc- 
ture and geometry of the nodes of A(p^) on the Fermi surface. 
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Order Parameter Models and Low Temperature Asymptotics 

Earlier experimentaSS and theoretical© II' @0'E1 transport studies 
indicated that the excitation gap in UPt3 has line nodes and probably polar 
point nodes. This is reinforced by the thermal conductivity data of Lussier 
et al.™ which is fit almost perfectly down to T ~ 50 mK by either the Ei g 
or the F>2u models, both of which have a line node in the aft-plane and polar 
point nodes. The differences in symmetry between the Ej s and E2« states are 
predicted to show up in the impurity dependences of the anisotropy ratio 
limT^o^c/ K 6 (see Fig. l).i However, as we show below the ground-state 
of the AE-model has cross nodes out of the basal plane located at polar 
angles, 9 = cos" _1 (±1/V3). The AE-model can also account for the thermal 
conductivity data along the c and b axes. The model also predicts a large 
anisotropy in the a6-plane. 



Tabic 1. Low temperature phases of several order parameter models. The first 
two entries are based on strong spin-orbit coupling and the symmetry group 
[Dsh] spin-orbit x T x £7(1). The third entry is based on no spin-orbit coupling 
and the -ESG-subclass SU{2) spin x [Dehlorbit x T x 17(1), and the last entry belongs 
to the mixed symmetry orbital subclass of the ESG, SO(3) sp i n - or bit x T x f7(l). 
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Generic basis functionsCJ for E-rep models, AB and ESG- models for 
the ground state of UPt3 are listed in Table 1. All of these models have 
both line nodes and point nodes - or "cross nodes" which are point nodes 
formed at the intersection of line nodes in the two different basis functions 
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of the accidental degeneracy models. Order parameter models like the polar 
state (Aiu), or combinations of the polar state with a nodeless Ai 9 state are 
excluded by existing transport data. 

Crystal symmetry alone does not determine the detailed structure of 
the order parameter or excitation gap as a function of position on the Fermi 
surface. In particular, the behavior of A( P y) in the vicinity of the nodes 
depends on material parameters. These 'nodal' parameters determine the 
rate at which the excitation gap opens near a linear point or line node (e.g. 
Ei 9 ground state), the curvature of a quadratic point node (e.g. F,2 U ground 
state) or the third derivative of a cubic point node (e.g. Bi n state). In 
principle these nodal parameters can be calculated by solving for the eigen- 
functions of the linearized gap equation, but this requires precise knowledge 
of the pairing interaction for momenta on the Fermi surface. Our approach is 
to enforce the symmetry of A( P j) through the positions of the nodes and to 
model the nodal structure with phenomenological parameters which we fix 
by comparison with experiment. Variational basis functions are constructed 
from the generic basis functions by multiplying by a function, Ty, which is 
invariant under all crystal group operations. This procedure allows us to 
introduce the Fermi surface anisotropy into the basis functions and indepen- 
dently model the nodal regions of the excitation spectrum for line and point 
nodes. Although the specific choice of functions Ty for a given representa- 
tion is not unique, once we fix the nodal parameters the low temperature 
thermal conductvity is determined. As we have shown previously,Bala the 
low-temperature response is dominated by the lower dimensional regions of 
the Fermi surface near the nodes of the order parameter; in these regions 
impurity induced Andreev states with energies e < 7 <C Ao form a 'metallic' 
band that determines the transport properties at k^T < 7. 

A key feature required to account for heat transport in UPt3 is the 
existence of gapless excitations with both vj-c 7^ and vyb 7^ 0. Machida et 
alx2 proposed an ESG-model based on no spin-orbit coupling and a 3D order 
parameter in spin space, corresponding to the three possible spin-triplet pair 
states, all of which have the same orbital pair wave function. They assume a 
ID orbital function with A2 U or Bi M symmetry.EJ The corresponding orbital 
basis functions are 

A A -( PJ ) = 



^bi u (p/;m,M3), (5) 

with A(9) « Aofi n 9 n near the point node, and A(9) w Ao/j,(j) near the line 
node at <j) = 0. The Bi u function has a cubic point and three axial line 
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nodes. The distinctive feature of this model is the large density of low- 
energy excitations near the poles. The A2 U function has a sixth-order point 
node and six axial line nodes. Among the ^4B-models, the A2 U ©Bi u order 
parameter, A AB (pf) = A A2u (pf) + iA Blu (p/), has the same nodal features 
as the Bi u model, i.e. a cubic point node and three axial line nodes. 

The density of impurity-induced Andreev states near the point and line 
nodes, and the heat current carried by these states, depends sensitively on 
the gap function near the point or line nodes. For a cubic point node at the 
poles, A(0 = 7r/6,0) »j ^AqO 3 , which leads to k c /T ~ /i 3 2 ^ 3 as T — > 0.0 
Thus, increasing ^3 tends to reduce the phase space for scattering among 
states near the cubic point node (see Fig. 2). This is the case for both Bi u 
and A2 U ©Bi„ models provided fj,s < 1. 




Fig. 2. The gap function |A(0)| sw /i3Ao|6*| 3 at the Fermi surface (FS) near a cubic 
point node. At low temperatures T < 7 <C Aq the excitations are confined to a 
lower dimensional region of the Fermi surface, i.e. < e < 7 near the node. The 
solid (dashed) line corresponds to a gap function with nodal parameter /X3 (yU 3 ). 
The dark shaded area is a measure of the increase in phase space for [13 < fi' a . 



For /i3 < 1 the low temperature limit of the anisotropy ratio of the 
thermal conductivity for the Bi u state is,c3 

^ MQ)M(o) 10 JL(&^) 1/S t (6) 



k c (T c ) I K b (T c ) /U 3 V 7 

where /x is the nodal parameter for the axial line nodes, i.e. A(<ft,0 = ir/2) ~ 
fiAo(j). For /i3 ~ 1 and 7/A0 <C 1 the ratio is typically much larger than 
unity. Based on Eq.(||) it appears possible to accomodate the experimental 
constraint, 7£ ex pt < 1, with a Bi„ order parameter by increasing ^3, i.e. 
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decreasing the phase space associated from the cubic point node.Il How- 
ever, Eq. (|6|) breaks down for /U3 > 1. This corresponds to shrinking and 
effectively removing the cubic point nodes. In this limit the axial line nodes 
determine the transport coefficients for both the c-axis and afc-plane cur- 
rents. We find that 7Z > 1, and in the limit /13 3> 1 the anisotropy ratio 
approaches unity (see also Table 2). 

The result 1Z > 1 also applies to the AB-model. This can be under- 
stood by noting that in the calculation of the thermal conductivity the 
order parameter enters as AA*. For the A2 M ffiBi u order parameter the 
product A A* = \A A2u (pf) + iA Blu (pf)\ 2 can always be written in the form 
I A Blu (p/ )| 2 (1 + J-(pf)), where T is an invariant function. This is precisely 
the variational form for the product for the basis function with Bi u sym- 
metry discussed above. Thus, we conclude that the k.2 U ®R\ u model cannot 
accomodate the experimental result 7£ e xpt < 1- The same arguments can 
be applied to the other relevant order parameter models of the ^45-class, 
e.g. Ai u ©Bi lt and k.\ g ®V>2 g - Thus, for a unitary order parameter belonging 
to the spin version of the i?5G-models (e.g. Bi„), and for the same parity 
subclasses of the ^45-models, e.g. A2 M ©Bi tt ), the anisotropy ratio is always 
greater than unity, which disagrees with the experimental results for 1Z£3 



Table 2. Asymptotic values of the thermal conductivity tensor n/T for T — > 0. 
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Machida et al. argue that UPt3 is described by a non-unitary 
spin-triplet ground state. Non-unitary ground states break time-reversal 
symmetry in spin space; S(pj) = iA(pj) x A*(p^) is proportional to 
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the spin polarization of triplet pairs with relative momenta pj. The 
ground state of the ESG- model of Ref. (13, 44) is a non-unitary state with 



A = A(p/)(a + ic)/2. The corresponding excitation spectrum has two 



branches: Ej = y£,p f + |A(p/)| 2 and Ej = \£, Pf \, where j (J.) corresponds 
to an excitation with spin parallel (anti-parallel) to S. Thus, the down-spin 
excitations are gapless over the whole Fermi surface and give rise to a large 
heat current at low temperatures compared to unitary or even parity su- 
perconducting states. The low temperature thermal conductivity is greater 
than half the normal-state value extrapolated below T c , i.e. k/T > | kjv/T. 
From Ref. (|26|) we deduce that the normal-state scattering rate extrapo- 
lated to T = is 1/t(0) « ~ X 1/t(T c ). This gives lim r ^ ^/T > [k n /T\ Tb , 
a result which is in striking contrast with the low temperature experimen- 
tal data for UPt3, which are roughly two orders of magnitude smaller, i.e. 
limT— »o K-i/T ~ 0.01[Kj/T]r c , with i = q.c. Thus, a non-unitary ground state 
of the form proposed by Machida et a/JiBEl is incompatible with the thermal 
conductivity data on UPt3.Q 

Now consider the orbital ESG-model proposed by Zhitomirsky and 
UedaHil They start from the five-dimensional 1 = 2 orbital representation of 
the full rotation group and assume that crystal-field terms break the full rota- 
tional symmetry and split the d-wave reprentation into three irreducible rep- 
resentations of D 6h ; r 2 [50(3)] => Ai 9 ©E l9 ©E 29 with T c Al9 > if 19 > if 29 . 
The onset of superconductivity is identified with the A% g component, and 
the second transition is identified with nucleation of an Ei 5 component. The 
AE-m.ode\ would in principle predict a third transition involving the E2 5 am- 
plitude, however, no such transition is observed in UPt3. Thus, T c 9 must 
be suppressed to very low temperatures by the development of the first two 
order parameter components. 

The onset of an Ei ff order parameter spontaneously breaks both time- 
reversal symmetry and rotational symmetry in the basal plane. The ground- 
state order parameter in the ^4i?-model has the form, A AE (pj) = A Al9 (~pf) + 
iA Elg (pf), with possible basis functions 



A Al9 ( P/ ) 



A^(T)i(3cos(0) 2 -i; 



FAigtPf, HA) , (7) 



A^( P/ ) = [Af(T)sin(2^sin(^)]^ l9 (p / ;^). (8) 

The Ai 9 basis function has tropical line nodes above and below the equator 
at 6 = cos _1 (±l/\/3)- The Ei g basis function has an axial line node run- 

1 Early transport and heat capacity experiments at higher temperatures on poorer quality 
crystals were interpreted in terms of a large residual density of states at the Fermi energy. 
Indeed this interpretation was the original argument for introducing the non-unitary spin- 
triplet pairing model for UPt3. 
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ning from the north pole to the south pole. Consequently the ground-state 
AE order parameter has four cross nodes located at the intersection of the 
tropical and axial line nodes. As can be seen from Table 2 the asymptotic 
values of the thermal conductivities of the ^4i£-model are nonuniversal^j 

An additional feature of the ^i^-model is that the low-temperature 
phase spontaneously breaks rotational symmetry in the aft-plane. The low- 
temperature phase of the Ginzburg-Landau (GL) functional for the AE- 
model is also degenerate with respect to rotations of the Ei 9 order parameter 
in the aft-plane. This accidental degeneracy is lifted by higher order terms in 
the free energy (sixth-order and higher) which differentiate hexagonal sym- 
metry from cylindrical symmetry. In order for the AE-model to account for 
the thermal conductivity data for heat flow along the c and 6 axes the ground- 
state must orient such that the axial line node of the Ei g order parameter 
is along the a-direction (A; x -direction). If the hexagonal anisotropy energy is 
small, J-hex ^ -^/|A| 2 > the heat current may orient the Ei g order parameter 
so that the current response is nearly isotropic in the aft-plane. However, for 
sufficiently small currents the anisotropy of the thermal conductivity tensor 
that results from the spontaneously broken rotational symmetry in the aft- 
plane should be observable. Table 2 summarizes the low temperature limits 
for the thermal conductitivities for the models discussed above. 



Numerical Results and Conclusions 

In order to obtain quantitative results for the thermal conductivity over 
the full temperature range below T c we must include inelastic 'electron- 
electron' scattering. We model the inelastic scattering rate in the super- 
conducting phase of UPt3 by extrapolating the normal-state scattering rate, 
l/ r ee = AT 2 , below T c . For T < 0.3T C the inelastic channel is unimportant 
even for the clean limit of UPt3 discussed here. The combined elastic and 
inelastic scattering rate in the normal-state is r(T) = Fq + AT 2 . From the 



measurements in Ref. 26 it follows that the UPt3 crystals that were studied 



2 

c > 



have nearly equal elastic and inelastic scattering rates at T c , i.e. Tq = AT C 
so that r(T) = ro(l + T 2 /T 2 ). Our procedure for computing and fitting 
the thermal conductivities is described in detail in Refs. |25| and S. Below 
T c we assume that the elastic scattering is in the unit arityl limit f Here we 
briefly outline the procedure and then discuss the numerical results for the 
Bi u , AB, and AE models. 

2 An extra complication in the analysis of the ^4-E-model is that the angular average of 
the order parameter is (generally) nonvanishing, (A Al9 (p/)} 7^ 0, which results in non- 
vanishing off-diagonal contributions to the scattering self-energy which must be calculated 
self-consistently. 
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We describe the analysis of the Bi u order parameter; a similar analysis 
is carried out for the other order parameter models. We first calculate the 
basal plane thermal conductivity, k&, and adjust the scattering rate To and 
the nodal parameters, \x at </> = 0, 9 = ir/2 and fi$ at <fi = ir/6,9 = 0, 
to obtain a good fit between theory and experiment. The next step is to 
compute Kf, and k c keeping Tq and u. fixed, and relaxing ^3 for a best fit. 
Figure 3 shows the nearly perfect fit to the basal plane thermal conductivity 
of Ref. (H) for the B lu model. 

For the variational function J-b 1u we have chosen the form 



•Fb\ u = 2/o(a<£, ae) {l + cos 2 (30) + a e cos 2 (6>)) , (9) 



where the prefactor yo is determined by maxA(pf) = Ao- The coefficients 
(a^ag) determine the nodal parameters, 



u. = 3(1 + a^) yoia^, a e ) , (10) 
/i 3 = (1 + a e ) yo(aj,,ae) ■ (11) 



Thus, we have two variational parameters (0^,0^) that enable us to inde- 
pendently adjust the nodal parameters /J, and ^3. 

We also examined the sensitivity of our numerical results to the form 
of the variational function. In particular, we calculated the thermal conduc- 
tivity for a piecewise continuous Bi u order parameter of the form A(0, 9) = 
Ajf lu g{<t>)h(9) with 



/iA o : < 4> < Y . ^ 



1 : 4>*<<t>< 



6 ' 



m - \ 1 : o,<e<%, (13) 



where = min(|, 1/^Ao), and 9* = min(|, l/( / U3Ao) 1 ^ 3 ), and the imposed 
B lu symmetry on the functions g(<p)h(9); i.e. g(^ —(/))= g(<ft), g(^ + <f)) = 
-g(4>), and h(ir -9) = h{9). 
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Fig. 3. (a) The normalized thermal conductivities for a Bi„ state in comparison 
with the experimental data on UPt 3 .E3 The best fit parameters are shown in 
the figure, (b) The same comparison for the A2 U © B\ u model, and (c) for a 
piecewise continuous Bi u state (see text). 
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With this piecewise Bi„ model we confirm our earlier conclusions, viz. 
that only the values of Ao, To, and {u^i} vary slightly, but the best fit low- 
temperature thermal response functions are unaffected. Furthermore, in the 
limit ^3 > 1 (ie. h(6) — > 1) we obtain TZ — ► 1 for the asymptotic anisotropy 
ratio (see Fig. 3c). In short, all the numerical results are in agreement 
with our estimates of the asymptotic values of the B\ u pair state. These 
considerations apply to the analysis of the AB and ESG-models as well. 

Figure 3 shows our 'best fit' calculations for the Bi„ (Fig. 3a) and 
AB (Fig. 3b) order parameter models. For both models we obtain excellent 
fits to the basal plane thermal conductivity, «&/T, over the entire temper- 
ature range. However, we are unable to fit the c-axis thermal conductivity 
with the Biu and AB order parameter models, which is consistent with our 
phase space arguments and asymptotic estimates. 

We obtain good low temperature fits for both the b and c axis thermal 
conductivity with the ^i^-model provided the Ei s component is oriented 
with the axial line node along the a-axis. The insert in Fig. 4 shows the 
strong anisotropy of the heat flow in the afe-plane for the AE-model, i.e. 
limj^o K a / Kb ~ 4, which is a direct consequence of the broken rotational 
symmetry in the basal plane. This anisotropy, which is specific to the the AE- 
model, has not been experimentally observed. It would serve as an important 
test of the ^.E-model and would differentiate the model from the 2D E-rep 
models. The nodal parameters, which determine the spectrum of low-lying 
excitations around the cross nodes, are defined as A A ^(9) « A$n A (9* - 9), 
with 0* = cos' 1 and A El f(0) « Af/x^. 

The quality of our fits for the thermal conductivity along different 
crystal axes is best summarized in terms of the anisotropy ratios 1Z(T) = 
[k c (T) I Kb(T)]/ (« C (T C ) / Kb(T c )] for the various order parameter models. Fig- 
ure 5 shows that the anisotropy ratio for the data of Ref. p6| ) is well de- 
scribed by the ground states for the Ei ff , E2 U and Ai 5 ©Ei 9 models. The AB- 
models, and the spin-triplet -ESG-models cannot account for the anisotropy 
of the thermal conductivity at low temperatures. 

Independent measurements of the thermal conductivity by Suderow 
et alE^ differ somewhat from those of Lussier et a/.,ll particularly in the 
anisotropy ratio 1Z. This difference is mainly due to differences in the c-axis 
thermal conductivity. The normalized basal plane conductivities, plotted vs. 
T/T c , are in remarkably good agreement, despite of being measured on dif- 
ferent samples. The differences in the data for heat flow along the c-axis for 
the different samples might be due to crystal imperfections along the c-axis, 
such as stacking faults, structural modulations or a mosaic structure. 
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Fig. 4. The normalized thermal conductivities for the Ai g ©Ei s model in com- 
parison with the experimental data on UPt3.il The inset shows the ratio of the 
conductivities Kt/K a . 




T/T 

c 

Fig. 5. The anisotropy ratio of the thermal conductivities for various 'best fit' 
pairing states. The lines for an Ei 3 , E2 U and Ai„ state are shown for comparison. 
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Conclusion 

Our analysis of the low temperature thermal conductivities for a single 
component order parameter of the ESG-class (e.g. Bi u ), and the accidental 
degeneracy order parameter models of the ^45-class (e.g. A2 U ©Bi u ), and 
the orbital subclass of the ESG-class (Ai s ®Ei g ), support earlier analyses 
suggesting that one of the 2D E-rep models is a promising model for the 
order parameter of UPt3. The AB- models and the spin-triplet version of the 
ESG-model are in substantial disagreement with the experimental data at 
low temperatures. Cubic point nodes or cross nodes at or near the poles 
of the Fermi surface lead to a large heat current along the c-axis, which is 
not reconcilable with the experimental data. The two-component orbital 
AE-m.ode\ can account for the c-axis and 6-axis thermal conductivities with 
an appropriate orientation of the Ei 5 order parameter. However, the AE- 
model spontaneously breaks rotational symmetry in the basal plane which 
so far has not been observed in any experiment to date. The basal plane 
anisotropy is large (lim^o i^a/^b ~ 4) and should easily be detected. 
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